Polarized quark distributions in nuclear matter 
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We compute the polarized quark distribution function of a bound nucleon. The Chiral Quark- 
Soliton model provides the quark and antiquark substructure of the nucleon embedded in nuclear 
matter. Nuclear effects cause significant modifications to the polarized distributions including an 
enhancement of the axial coupling constant. 
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Polarized lepton-nucleus scattering experiments are an 
important tool in hadronic physics. For example, in or- 
der to study the spin structure function of the neutron, 
one must use nuclear targets. It is already well known 
that there are significant differences betweeen free and 
bound nuclcons in the unpolarized case; the famous Eu- 
ropean Muon Collaboration (EMC) effect pj is the prime 
example. It is reasonable to assume that nuclear effects 
could appear in polarized quark distributions. Our pur- 
pose here is to calculate the analogous modification to 
the nucleon spin structure function g^' n \x,Q 2 ): a 'po- 
larized EMC effect'. 

The first discussion of nuclear effects in the polarized 
quark distributions is in Ref. in the context of dy- 
namical rescaling. A more recent calculation Q pre- 
dicts dramatic effects for the bound nucleon spin struc- 
ture function. We have shown 0, 0] that sea quarks, 
introduced at the model scale, can have important con- 
sequences for modifications in the nuclear medium. We 
will use our previous work 4| as a basis for the results 
presented here, and provide a mechanism for the mod- 
ification within the Chiral Quark-Soliton (CQS) model 
IH IE 111 113 • This relativistic mean field approxima- 
tion to baryons has many desirable qualities such as the 
inclusion of antiquarks (which is deeply linked to satis- 
fying sum rules and the positivity of Generalized Parton 
Distributions), and a basis in QCD ||. We have previ- 
ously shown how the model describes nuclear saturation 
properties, reproduces the EMC effect, and satisfies the 
bounds on unpolarized nuclear antiquark enhancement 
provided by Drell-Yan experiments |4|. Therefore, we 
expect the CQS model to produce a reasonable result for 
the polarized distributions. 

The CQS model Lagrangian with (anti)quark fields 
V*,^, and profile function 0(r) is 



C = tp(ip - Me 



(i) 



where 0(r — ► oo) = and 0(0) = — n to produce a 
soliton with unit winding number. The quark spectrum 
consists of a single bound state and a filled negative en- 
ergy Dirac continuum; the vacuum is the filled negative 
continuum with = 0. In both the free nucleon and vac- 
uum sectors the positive continua are unoccupied. The 
wave functions in this spectrum provide the input for the 
quark and antiquark distributions used to calculate the 



nucleon structure function. 

We work to leading order in the number of colors 
(N c = 3), with Nf = 2, and in the chiral limit. While 
the former characterizes the primary source of theoretical 
error, one could systematically expand in Nq to calculate 
corrections. We also expect that since the nucleon size 
is stable in the limit Nc — > oo, the quark wavefunctions, 
our primary fo cus, should be within a few percent of their 
Nc = 3 value [Tl|. We take the constituent quark mass 
to be M — 0.42 GeV, which reproduces, for example, 
the A-A mass splitting at higher order in the Nc expan- 
sion, and other observables We ignore contributions 
from the structure functions of pion quanta, which in 
this model propagate through constituent quark loops; 
they are suppressed by factors of 0(1/Nc), and are not 
treated at leading order. 

The theory contains divergences that must be regu- 
lated. We use a single Pauli-Villars subtraction as in 
Ref. 01 because we follow that work to calculate the 
quark distribution functions. The Pauli-Villars mass is 
determined by reproducing the measured value of the 
pion decay constant, f„ = 0.093 GeV, with the rele- 
vant divergent loop integral regularized using Mpy — 
0.58 GeV. This regularization also preserves the com- 
pleteness of the quark states [Pj • 

The results for binding and saturation of nuclear mat- 
ter have been published elsewhere but we provide a 
brief review for completeness. The nucleon mass is given 
by a sum of the energy of a single valence level (E v ), 
and the regulated energy of the soliton (Eq equal to the 
energy in the negative Dirac continuum with the energy 
in the vacuum subtracted) 



M N = N c E v + E e (M)~ 
E e (M) = Nc^2 E n- E n ] 



M 2 
M 2 ^ 



Ee(M PV ) (2) 
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The field equation for the profile function is 
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(3) 



(4) 



where p% and p^ s are the quark scalar and pseudoscalar 
densities, respectively. The dependence of nucleon prop- 
erties on the nuclear medium has been incorporated in 
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the model by simply letting the quark scalar density in 
the field equation Q contain a constant, but Fermi mo- 
mentum kp dependent, contribution, P^(k F ), equal to 
the convolution of the nuclear scalar density with the 
nucleon quark density 



P?(k F )= / d 3 r'p»(r')pl(r-r>) 



(5) 



arising from other nucleons present in symmetric nuclear 
matter. This models a scalar interaction via the exchange 
of multiple pairs of pions between nucleons, and the pa- 
rameter g s is varied to obtain nuclear saturation. 

The nucleon scalar density is determined by solving 
the nuclear self-consistency equation 



= 4 



d 3 k 



M N (p : 



(2tt) 3 y/k* + M N (p?)* 



(6) 



The dependence of the nucleon mass, and any other prop- 
erties calculable in the model, on the Fermi momentum 
k F enters through Eq. 10. Thus there are two coupled 
self-consistency equations: one for the profile, Eq. J3J, 
and one for the density, Eq. ©. These are iterated un- 
til the change in the nucleon mass Eq. J5J is as small 
as desired for each value of the Fermi momentum. We 
use the Kahana-Ripka (KR) basis 01 to evaluate the 
energy eigenvalues and wave functions used as input for 
the densities, nucleon mass, and quark distributions. 

We introduce a phenomenological vector meson (with 
mass fixed at m v = 0.77 GeV and coupling g v ) [l4[ ex- 
changed between nucleons, but not quarks in the same 
nucleon (i.e. we ignore the spatial dependence of the vec- 
tor field in the vicinity of a nucleon, treating only the 



nuclear mean field), 
vector density 



The vector meson couples to the 



P. 



N 



(k F ) = J dVpf(r>*(r 



9k 3 
' 3tt 2 ' 



(7) 



This mechanism is a proxy for uncalculated soliton- 
soliton interactions used to obtain the necessary short 
distance repulsion which stabilizes the nucleus. 

The polarized quark distribution for flavor i is defined 
by the difference between the quark distributions with 
spin parallel (j) and antiparallcl (J.) to the nucleon 



A qi (x,Q 2 ) = qj(x,Q 2 )-qi(x,Q 2 ). 



(8) 



The polarized antiquark distribution is defined analo- 
— t — i 

gously using q\ , and q\ . The isovector polarized distri- 
bution Aq^ T=v > '(a;) — Au(x) — Ad(x) is the leading order 
term in Nq, with the isoscalar polarized quark distribu- 
tion Aq( T=0 ) (x) = Au(x) + Ad(x) smaller by a factor 
~ 1/Nc and set to zero. This follows from the fact that 
the isoscalar combination is normalized to the spin of the 
nucleon, which is O(Nq), while the isovector combination 
is normalized to the axial coupling, which is 0(Nq) [l5| . 
Therefore, at the model scale Mp v ~ 0.34 GeV 2 , we see 
that a large portion of the spin is carried by the orbital 
motion of the constituent quarks in the valence level and 
the sea . We will therefore suppress the isospin super- 
script in the following. The distributions are calculated 
using the KR basis at kp — and kp = 1.38 fm _1 (see 
Refs. 0,O|) almost exactly as in Ref. ^3 where the quark 
distribution is given by the matrix element 



Ag(x) = -\(2T 3 )N c M N Y / (^n\r 3 l5 (l + j 1 3 )S(E n + P 3 - xM N )\iP ri 



(9) 



r 



with the regulated sum taken over occupied states. The 
eigenvalues E n are determined from diagonalizing the 
Hamiltonian, derived from the Lagrangian l|T[l. in the 
KR basis. These are also the eigenvalues that enter into 
Eq. J2J) for the mass. The momentum sum rule (for 
the unpolarized distribution) is automatically satisfied 
as long as Eq. © defines the mass in the unpolarized 
analog of Eq. @, and the same eigenvalues are used in 
both equations |12l |. 

The vector and scalar interactions with the other nu- 
cleons in the nucleus at the (low) model scale are implic- 
itly included in the energy eigenvalues in Eq. ©. In the 
'handbag' diagram for deep inelastic scattering in the 
parton model language, the quarks in the intermediate 
state are treated as non-interacting, so how do the nu- 
clear interactions modify the parton distributions? The 



key point is that all three quarks in the intermediate state 
undergo evolution in QCD from the same starting scale 
to the scale of the deep inelastic scattering, and it is this 
scale that would appear in the Wilson coefficients in the 
language of the operator product expansion (OPE) pic- 
ture of deep inelastic scattering. The model for medium 
modifications presented here represents different bound- 
ary conditions on the Wilson coefficients or the parton 
distributions for free and bound nucleons at the model 
scale that maintains consistency with the parton model 
and OPE pictures of deep inelastic scattering. It is worth 
noting here that these two pictures are already consistent 
in the free nucleon case since the parton model hypothesis 
that the quark transverse momenta do not grow with Q 2 
is satisfied jl'J , and our model for medium modifications 
does not damage this equivalence. 
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The antiquark distribution is given by Aq(x) — 
Aq(—x) where the sum is over unoccupied states. The 
use of a finite basis causes the distributions to be discon- 
tinuous. These distributions are smooth functions of x in 
the limit of infinite momentum cutoff and box size, but 
numerical calculations are made at finite values and leave 
some residual roughness. This is overcome in Ref. [T^| by 
introducing a smoothing function. We deviate from their 
procedure, and do not smooth the results; instead we find 
that per forming the one-loop perturbative QCD evolu- 
tion provides sufficient, but not complete, smooth- 
ing. Some residual fluctuations due to the finite basis 
remain visible in our results, and the size of these fluctu- 
ations serve as a guide to the size of the error introduced 
by the method. 

These distributions are used as input at the model 
scale of Q 2 = M PV ~ 0.34 GeV 2 for evolution to 
Q 2 = 10 GeV 2 . The polarized structure function to lead- 
ing order in Nq is given by 



Hi 



±18** 



NS 



0(N° C ) 



(10) 



(11) 



Aq NS = -(Am- Ad + Au- Ad) +0{N°). (12) 
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The ratio function is defined to be 



Ri{x,Q 2 ) 
g[ plA \x,Q 2 ,k F ) 



g^ A \x,Q 2 ,k F ) 



Ag[ p \x,Q2,k f 



0) 



(13) 



/ ^f(y)g[ p) (x/ yi Q 2 ,k F ). 
j x y 



The nucleon momentum distribution f(jj) in light polar- 
ized nuclei has been calculated in Ref. [I3 • Here, the nu- 
cleon momentum distribution is assumed to be the same 
as the unpolarized case, £ls the effects of the spin-orbit 
force will tend to average out in nuclear matter. We can 
also justify this approximation in nuclear matter because 
the zero pressure condition P + = P~ for a nucleus with 
momentum P in the rest frame, which implies the light- 
cone version of the Hugenholtz-van Hove theorem jig , 
is still true. Therefore, one expects a distribution f(y) 
that is peaked at y ~ 1, like those in Ref. 17]. This peak 
location is the dominant effect on the ratio Eq. (|13fl ; the 
remaining details of the function f(y) have only a small 
effect. Following a light-cone approach valid for any mean 
field theory of nuclear matter for which the density and 
binding energy per nucleon are the only input parameters 
[l8| one obtains 



f(y) 



4A|/ 



(l+A F -y)0(y-l + A F ) [A 2 , - (1 - y) 2 } 



_ _ (14) 

where A F = k F /M N and M N = M N (0) - 15.75 MeV. 

We show the ratio Eq. (|13fl in Fig. fusing a 'valence'- 
like distribution, as well as for the full distribution. The 
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FIG. 1: The ratio Eq. JTHJ at scale Q 2 = 10 GeV 2 for nuclear 
matter. The heavy line is the full calculation for nuclear mat- 
ter. The light line is the effect calulated using only medium 
modifications to the 'valence' energy level as decribed in the 
text. 



latter includes all medium modifications, while the for- 
mer distribution uses the medium modified energy level 
eigenstate, but the same free nucleon sea quark distribu- 
tion for both the free and bound nucleon. This was done 
in order to compare our results with the model in Ref. , 
which only has valence quarks at the model scale. The 
single energy level actually has a contribution to the po- 
larized antiquark distribution, so it alone cannot be con- 
sidered a true valence spin structure function. However, 
this contribution is small, so we effectively reproduce the 
result of a valence quark model, especially in the region 
x>0.3. 

In Fig.^ one can see that there is a large depletion for 
0.3 < x < 0.7 in the polarized 'valence' quark distribu- 
tion. This produces a large depletion in the isovector ax- 
ial coupling g^ of 17.8%. This large effect is comparable 
to that of the calculation in Ref. [3j which only includes 
valence quarks at the model scale. This valence effect is 
mitigated by a large enhancement in the sea quark con- 
tribution, so that the full polarized distribution has only 
a moderate depletion in the region 0.3 < x < 0.7 of the 
same size as the EMC effect in unpolarized nuclear struc- 
ture functions. There is a large enhancement for x < 0.3 
due to the sea quarks. This large enhancement is very dif- 
ferent from the small effect calculated in the unpolarized 
case 0, and seen in unpolarized Drell-Yan experiments 
[l9j . This would suggest that one might see a significant 
enhancement in a polarized Drell-Yan experiment, even 
after including shadowing corrections (which we address 
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later). The larger sensitivity to the lower components of 
the wave functions is the primary source for the greater 
sea quark enhancement in the polarized case, in contrast 
to the unpolarized case. 

The axial coupling g A ' is enhanced by 9.8% in the 
nuclear medium. This is in accord with an earlier finding 
of a ~ 25% enhancement for gA in a different soliton 
model by Birse [23 . There, the effect is also seen as 
a competition between enhancement and depletion. In 
order to address the medium modification of the Bjorken 
sum rule |2lll22l| 



lim 



dxg¥\x,Q*)- 9 r(x,Q*) = Y (15) 



» 



9A 



as an integral of the experimentally observed nuclear dis- 
tribution, one must account for the effects of shadowing. 
This occurs when the virtual photon striking the nucleus 
fluctuates into a quark-antiquark pair over a distance 
~ \/2Mfqx exceeding the inter-nucleon separation. This 
causes a depletion in the structure function for x < 0.1 
and is relatively well understood |23l |24j, |25j . Shadowing 
in the polarized case is expected to be larger than in the 
unpolarized case by roughly a factor of 2 simply from the 
combinatorics of multiple scattering (see e.g. Ref. 26]). 

The enhancement at x ^ 0.1 — 0.2 in Fig. ^is com- 
parable to that seen by Guzey and Strikman |2fj; they 
assume that the combined effects of shadowing, enhance- 
ment, and target polarization lead to the empirical value 
of the nuclear Bjorken sum rule for 3 He and Li. Shadow- 
ing effects become large for x < 0.05, but we ignore them 
as well as target polarization; such precision is not nec- 
essary for our relatively qualitative analysis. One needs 
~ 10 times the shadowing observed in the unpolarized 
case for Lead in order to counter the enhancement at 
x ~ 0.1 — 0.2, and give the same value for the Bjorken 
sum rule (|T5)l in matter and free space. This assumes 
that shadowing is the only significant effect neglected at 
small x in our calculation of the unpolarized quark dis- 
tribution Q. 

We also present, in Fig. [3 the results for the spin asym- 
metry 



E 



>,Q 2 )+ft(z,Q 2 )] 



(16) 



The nuclear asymmetry A^ A ^ is defined by replacing 
the polarized and unpolarized quark distributions, repre- 
sented generically as q, with 



^-f(y)q^(x/y,Q 2 ,k F ). 



(17) 
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FIG. 2: The asymmetry A[ plA) Eq. d) at scale Q 2 = 
10 GeV 2 . The heavy line is for nuclear matter. The dashed 
line is for the free proton. The data are for the free proton 
from SLAC (filled) for Q 2 ~ 1 - 40 GeV 2 and HERMES 
[H (empty) for Q 2 ~ 1-20 GeV 2 The free curve falls slightly 
below the data due to the lower value of gA calculated in the 
large Nc limit. 



The central mechanism to explain the EMC effect is 
that the nuclear medium provides an attractive scalar 
interaction that modifies the nucleon wave function. We 
see this again in the polarized case. This is also the 
dominant mechanism in the model of Cloet et al Q , and 
the soliton model of Birse [2(j ■ 

The present model provides a intuitive, qualitative 
treatment that maintains consistency with all of the free 
nucleon properties calculated by others H,@. It provides 
reasonable description of nuclear saturation properties, 
reproduces the EMC effect, and satisfies the constraints 
on the nuclear sea obtained from Drell-Yan experiments 
with only two parameters for the nuclear physics (g s 
and g v ) fixed by the binding energy and density of nu- 
clear matter. Therefore, we expect the results presented 
here to manifest themselves in future experiments with 
polarized nuclei. Our conclusions differ from those in 
Ref. 0; the main difference is the role of sea quarks at 
the model scale. Therefore, we also expect future exper- 
iments would help determine the role of sea quarks in 
nuclei. 



We find that for the free case, the calculation falls slightly 
below the data due to the smaller value of gA in the large 
Nc limit, and that the size of the medium modification 
is of the same order as the experimental error for the free 
proton [27LI2H. 
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